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The first direct detection limits on dark matter in the MeV to GeV mass range are presented, 
using XENON10 data. Such light dark matter can scatter with electrons, causing ionization of atoms 
in a detector target material and leading to single- or few-electron events. We use 15 kg-days of data 
acquired in 2006 to set limits on the dark-matter-electron scattering cross section. The strongest 
bound is obtained at 100 MeV where a e < 3 x 10 -38 cm 2 at 90% CL, while dark matter masses 
between 20 MeV and 1 GeV are bounded by o e < 10~ 37 cm 2 at 90% CL. This analysis provides a 
first proof-of-principle that direct detection experiments can be sensitive to dark matter candidates 
with masses well below the GeV scale. 



INTRODUCTION. Most current dark matter (DM) 
direct detection experiments focus on detecting a Weakly 
Interacting Massive Particle (WIMP) with a mass of 
1 — 1000 GeV. There are two main reasons for this focus. 
Theoretically, a WIMP in this mass range can naturally 
have the correct thermal relic abundance [T] . Experimen- 
tally, a WIMP-nucleus scattering event is likely to pro- 
duce detectable quanta (phonons, scintillation photons, 
ionization or some combination of these). DM candidates 
with mass < 1 GeV typically cannot produce nuclear re- 
coil signals above detector thresholds, and have therefore 
been largely ignored. 

It is straightforward, however, to theoretically con- 
struct well-motivated, viable, and natural DM candidates 
with sub-GeV masses (e.g. [SIS])- Given the current lack 
of firm experimental evidence for WIMPs in any mass 
range, it is important to search for other theoretically mo- 
tivated DM candidates. As was recently proposed in [2], 
sub-GeV DM can lead to observable signals if it scatters 
with atomic electrons, as opposed to nuclei. This scat- 
tering can ionize atoms in a target material, resulting in 
single-electron signals. As discussed below, few-electron 
signals may result if the primary ionized electron or de- 
excitation photons lead to further ionization. 

Dual-phase liquid xenon detectors have demonstrated 
sensitivity to such small ionization signals [7HS]- In this 
letter, we present the first direct detection limits on 
MeV-GeV-mass DM, using 15 kg-days of exposure of the 
XENON10 experiment obtained with a single-electron 
trigger threshold [10]. We consider the observed rate of 
one-, two-, and three-electron events. The origin of these 
events is unclear, and they are likely to result from back- 
ground processes. The data nevertheless allow robust 
limits to be set for DM as light as a few MeV. 

DATA SAMPLE. The XENON10 Collaboration has 
reported results from a 12.5 live-day search for scatter- 
ing of low-mass (few-GeV rather than sub-GeV) WIMPs 
with xenon nuclei |10j . Particle interactions in the liq- 



uid xenon target can produce both ions (Xe + ) and ex- 
cited atoms (Xe*). A fraction of the ions recombine to 
form other Xe*, whose de-excitation process produces 
7eV scintillation photons. Electrons that escape recom- 
bination are accelerated away from the interaction site by 
an electric field, and extracted from the liquid to the gas 
with an efficiency that is essentially unity [TTJ [T^] ■ Under 
the influence of a high electric field in the gaseous xenon 
(~10kV/cm), each extracted electron produces 0(100) 
scintillation photons [13]. The detector's array of pho- 
tomultiplier tubes measures an average of 27 of these 
photoelectrons per extracted electron. 

The search for few-GeV dark matter reported in [TU] 
imposed a software ionization threshold of 5 electrons due 
to uncertainties in the ionization yield from very low- 
energy nuclear recoils. However, the XENON10 hard- 
ware trigger (described in [9]) was sensitive to single 
electrons. Among liquid xenon targets, this is the lowest 
trigger threshold in all reported dark matter search data. 
However, the precise trigger efficiency for this data sam- 
ple was not reported. To better understand the trigger 
efficiency, we have performed a detailed Monte Carlo sim- 
ulation of the trigger response to single to few-electron 
events in the XENON10 detector, based on the infor- 
mation given in Sec. 2.8 of [Hj. This simulation allows 
the hardware trigger efficiency to be calculated, and the 
result is shown in Fig. [I] ( bottom) . 

The performance and accuracy of our trigger efficiency 
simulation has been verified by comparing its prediction 
to the observed trigger roll-off of the calibration spec- 
trum reported in |15j , whose conditions differed from the 
present data only in the hardware threshold set point. 
The good agreement in this known case confirms the va- 
lidity of the simulation, which is then left with a single 
free parameter: the hardware threshold set point. We 
constrain this threshold by noting that the trigger effi- 
ciency curve must "turn-on" at, or prior to, the first non- 
zero bin in the measured spectrum of triggering events, 
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FIG. 1: Top: The spectrum of XENON10 dark matter 
search data, corrected for trigger efficiency. Blue boxes in- 
dicate statistical uncertainty, while green boxes indicate the 
systematic uncertainty arising from the the trigger efficiency. 
The efficiency curve crosses 5% within the orange-hatched ver- 
tical band. The thick gray curve is the best-fit triple Gaussian 
function. Thin solid red curves indicated the best-fit individ- 
ual components. Dashed lines indicate curves allowed at the 
90% upper limit for each component. Small open squares in- 
dicate the raw spectrum (uncorrected for trigger efficiency) 
from [10]. Arrows indicate 1-cr upper limits on the number 
of events for bins with no events. Bottom: The trigger effi- 
ciency as determined by Monte Carlo simulation, whose range 
is chosen such that the efficiency curve crosses 5% at, or be- 
fore, the first non-zero bin in the blue histogram. 

shown in Fig. 2 of [10] ■ In this context, we define the 
turn-on point as the location where the efficiency curve 
crosses 5%, which is indicated by the orange-hatched ver- 
tical band in Fig. [TJ If the efficiency were to turn on at a 
higher point, the peak of the single-electron distribution 
would be shifted to values much lower than that of the 
known detector response to these events, demonstrated 
by Fig. 2 (top) of [TO]. 

The measured spectrum of triggering ionization events, 
which we analyze for a signal, is given in Fig. 2 (top) 
of [TDJ- We reproduce this spectrum in Fig. [TJ (top), cor- 
rected for the trigger efficiency. Wide (blue) bars rep- 
resent statistical uncertainty, while the narrow (green) 
bars indicate the systematic uncertainty introduced by 
the range of allowed trigger efficiencies. This spectrum 
is fit by a triple Gaussian function with five free pa- 
rameters: the heights, Hi, of the three components and 
the mean and width of the first component (fix, ci). 
The means, [ii, and widths, o~i, are constrained to fol- 
low the relations [ii — and er,; = a\\fi, respectively, 
where i = 1,2,3 identifies the Gaussian component. In- 
dividual marginal posterior probability distributions are 
obtained for the event rates of the three components, 



ri = H i a i \'2iT/eSAx, where e = 0.92 is the overall cut 
efficiency reported in [ID] . S = 15 kg-days is the exposure, 
and Ax = 0.1 electrons is the histogram bin width. From 
these, upper limits are extracted taking the measured 
spectrum to be due entirely to signal (i.e. no background 
subtraction). The result of the fit, including statisti- 
cal and systematic uncertainties, gives 90% upper con- 
fidence bounds of r\ < 23.4, r 2 < 4.23, and r 3 < 0.90 
cts kg" 1 day" 1 . 

DIRECT DETECTION RATES. We assume that 
DM particles scatter through direct interactions with 
atomic electrons. If the DM-electron interaction is in- 
dependent of the momentum transfer, q, then it is com- 
pletely parametrized by the elastic cross section, <r e , of 
DM scattering with a free electron. For g-dependent in- 
teractions, we define a cross section a e by fixing q = am e 
in the matrix element [2] ■ The q dependence of the matrix 
element is then described by a DM form- factor, Fdm(9); 
for example, if the interaction proceeds through a mass- 
less vector mediator then Fdm = (am e /q) 2 . 

A large fraction of the kinetic energy carried by a DM 
particle, Fdm = mDMf 2 /^ — 10 eV(mDM/20 MeV), can 
be transferred to a primary ionized electron. We treat the 
target electrons as single-particle states bound in isolated 
xenon atoms, using the numerical RHF bound wavefunc- 
tions tabulated in [IB]. The electron recoils with energy 
F or , with a differential ionization rate [2] 



dRi, 



d\nE r 



Nn 



Pbm 
m DM 



E 

nl 



d(a 



ill 



d\nE r 



(1) 



where Nt is the number of target atoms, pdm 
0.1 GeV cm -3 is the local DM density, and the velocity- 
averaged differential ionization cross section for electrons 
in the (n, I) shell is given by 
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E et )/q + q/2m BM , and T)(v mi ^ 
has its usual meaning (^9(v— UmiJ)- We assume a stan- 
dard Maxwell-Boltzmann velocity distribution with cir- 
cular velocity vo = 220 km s" 1 and a hard cutoff at 
Vcsc = 544 km s" 1 P2]. 

With full shells, the form-factor for ionization of an 
electron in the (n, I) shell, escaping with momentum k' = 
\/2m e E CT after receiving a momentum transfer q, can be 
written as 
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where [■ ■ • ] is the Wigner 3-j symbol and Jl is a spher- 
ical Bessel function. The radial wavefunctions Rk'i'ir) 
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of outgoing electrons are found by numerically solving 
the radial Schrodinger equation with a central potential 
Zes(r)/r. Z c ff(r) is determined from the initial electron 
wavefunction, assuming it to be a bound state of the same 
central potential. We evaluate the form-factors numeri- 
cally, cutting off the sum at large I', L once it converges. 
Only the ionization rates of the 3 outermost shells (5p, 
5s, and 4d, with binding energies of 12.4, 25.7, and 75.6 
eV, respectively) are found to be relevant. 

The energy transferred to the primary ionized electron 
by the initial scattering process is ultimately distributed 
into a number of (observable) electrons, n e , (unobserved) 
scintillation photons, n 7 , and heat. To calculate n e , we 
use a probabilistic model based on a combined theoreti- 
cal and empirical understanding of the electron yield of 
higher-energy electronic recoils. Absorption of the pri- 
mary electron energy creates a number of ions, Ni, and 
a number of excited atoms, N cx , whose initial ratio is 
determined to be N ex /Ni « 0.2 over a wide range of ener- 
gies above a keV [HI [19]. Electron-ion recombination ap- 
pears well-described by a modified Thomas-Imel recombi- 
nation model [501 [5T] i which suggests that the fraction of 
ions that recombine, f R , is essentially zero at low energy, 
resulting in n e = Ni and n 7 = iV ex . The fraction, f e , 
of initial quanta observed as electrons is therefore given 
by f e = (1 - /b)(1 + AW^r 1 « 0.83 [21, The total 
number of quanta, n, is observed to behave, at higher 
energy, as n = E eI /W, where E cr is the outgoing energy 
of the initial scattered electron and W = 13.8 eV is the 
average energy required to create a single quanta [23] . 
As with fx and N ex /Ni, W is only well measured at en- 
ergies higher than those of interest to us, and thus adds 
to the theoretical uncertainty in the predicted rates. We 
use N cx /N t = 0.2, f R = and W = 13.8 eV to give 
central limits, and to illustrate the uncertainty we scan 
over the ranges < f R < 0.2, 0.1 < N cx /N, < 0.3, 
and 12.4 < W < 16 eV. The chosen ranges for W and 
N ex /Ni are reasonable considering the available data 
[9] HH1 HH H2]- The chosen range for f R is conserva- 
tive considering the fit of the Thomas-Imel model to low- 
energy electron-recoil data [20] . 

We extend this model to DM-induced ionization as fol- 
lows. We calculate the differential single-electron ion- 
ization rate following Eqs. ([l]]3]). We assume the scat- 
tering of this primary electron creates a further rS 1 ^ = 
Floor(E er /W) quanta. In addition, for ionization of the 
next-to-outer 5s and 4d shells, we assume that the pho- 
ton associated with the de-excitation of the 5p-shell elec- 
tron, with energy 13.3 or 63.1 eV, can photoionize, cre- 
ating another rfi 2 ^ = (1) or 4 quanta, respectively, for 
W > 13.3 eV (< 13.3 eV). The total number of detected 
electrons is thus n e = n' e + n", where n' e represents the 
primary electron and is thus or 1 with probability f R 
or (1 — f R ), respectively, and n" follows a binomial dis- 
tribution with n' 1 ) + n^ 2 ) trials and success probability 
f e . This procedure is intended to reasonably approxi- 
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FIG. 2: Top: Expected signal rates for 1-, 2-, and 3-electron 
events for a DM candidate with a e — 10 -36 cm 2 and Fdm = 1. 
Widths indicate theoretical uncertainty (see text). Bottom: 
90% CL limit on the DM-electron scattering cross section 
cr e (black line). Here the interaction is assumed to be in- 
dependent of momentum transfer (Fdm = 1)- The dashed 
lines show the individual limits set by the number of events 
in which 1, 2, or 3 electrons were observed in the XENON10 
data set, with gray bands indicating the theoretical uncer- 
tainty. The light green region indicates the previously allowed 
parameter space for DM coupled through a massive hidden 
photon (taken from [2]). 

mate the detailed microscopic scattering processes, but 
presents another 0(1) source of theoretical uncertainty. 
The 1-, 2-, and 3-electron rates as a function of DM mass 
for a fixed cross section and -Fdm = 1 are shown in Fig. [2] 
(top). The width of the bands arises from scanning over 
/h> N cx /Ni and W, as described above, and illustrates 
the theoretical uncertainty. 

RESULTS. Fig. [2] (bottom) shows the exclusion limit in 
the mcM-fe plane based on the upper limits for 1-, 2-, 
and 3-electrons rates in the XENON10 data set (dashed 
lines), and the central limit (black line), corresponding 
to the best limit at each mass. The gray bands show the 
theoretical uncertainty, as described above. This bound 
applies to DM candidates whose non-relativistic inter- 
action with electrons is momentum-transfer independent 
(Fdm = !)■ For DM masses larger than ~15MeV, the 
bound is dominated by events with 2 or 3 electrons, due 
to the small number of such events observed in the data 
set. For smaller masses, the energy available is insuffi- 
cient to ionize multiple electrons, and the bound is set 
by the number of single-electron events. The light green 
shaded region shows the parameter space spanned by 
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FIG. 3: 90% CL exclusion limits on a e for candidates with 
a DM form-factor -Fdm(?) of (am e /q) (red/lower line), cor- 
responding to DM with an electric dipole moment, and 
(am.e/q) 2 (blue/upper line), corresponding to DM scattering 
through a very light mediator. Dashed lines and bands are as 
in Fig. [2] The pale blue region shows the previously allowed 
parameter space for DM coupled through a very light hidden 
photon (Fdm = (am e /q) 2 ), with the gray strip indicating the 
"freeze-in" region (taken from [5]). 

models in which the DM candidate is a fermion coupled 
to the visible sector through a kinetically mixed "hid- 
den photon" with O(MeV-GeV) mass, and satisfying all 
previously known constraints (from [2]; see also [3] [25]). 

Fig. [3] shows the exclusion limits in the muM-Pe plane 
for DM candidates whose interaction with electrons is 
enhanced at small momentum-transfers by a DM form- 
factor, Fum ■ The red (lower) curves correspond to 
Fdm — (& m e/q), or DM scattering through an elec- 
tric dipole moment, and the blue (upper) curves to 
Fdm — {am e /q) 2 , or DM scattering though a very light 
(<§; keV) scalar or vector mediator. Bounds set by 1-, 
2-, and 3-electron rates are shown by dashed lines, and 
the central limits by dark lines. The bands illustrate 
the theoretical uncertainty. Both form-factors suppress 
the relative rate of events with larger energy deposition, 
and so reduce the fraction (and hence the importance) of 
events containing multiple electrons. The pale blue re- 
gion shows the parameter space for DM coupled through 
a very light hidden-photon mediator, and satisfying all 
previously known constraints, with the gray strip show- 
ing where the correct abundance is achieved through 
"freeze-in" (from j2]). These regions should be compared 
to the blue exclusion curve. 

DISCUSSION. The results above demonstrate, for the 
first time, the ability of direct detection experiments to 
probe DM masses far below a GeV. It is encouraging that 
with only 15 kg-days of data, and no attempt to control 
single-electron backgrounds, the XENON10 experiment 
places meaningful bounds down to masses of a few MeV. 

It should be emphasized that this analysis lacks the 
ability to distinguish signal from background. One 



promising method is the expected annual modulation of 
the signal. As discussed in [2], additional discrimination 
may be possible via the collection of individual photons, 
phonons [23], or ions, although at present such technolo- 
gies have yet to be established. 

Independently, this type of search could be signifi- 
cantly improved with a better understanding of few- 
electron backgrounds. A quantitative background es- 
timate was not made in [10] . making background sub- 
traction impossible. Single-electron ionization signals 
have been studied, and potential causes discussed, by 
XENON10 [2], ZEPLIN-II [7j, and ZEPLIN-III gj. Pos- 
sible sources include photo-dissociation of negatively 
charged impurities, spontaneous emission of electrons 
that have become trapped in the potential barrier at the 
liquid-gas interface, and field emission in the region of 
the cathode. The former two processes would not be 
expected to produce true two- or three-electron events, 
although single electron events may overlap in time, giv- 
ing the appearance of an isolated, double-electron event. 
With a dedicated study, these backgrounds could be 
quantitatively estimated and reduced. 

With larger targets and longer exposure times, ongo- 
ing and upcoming direct detection experiments such as 
XENON100, XENON1T, LUX, and CDMS, should be 
able to improve on the sensitivity reported here. Such 
improvements may require optimizations of the trigger- 
ing thresholds, and will strongly benefit from additional 
studies of the backgrounds. 
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